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Abstract. The large difference between the Planck scale and the electroweak scale, known as the hierarchy 
problem, has been addressed in some models through the existence of extra spatial dimensions. A search for 
evidence of extra spatial dimensions has been performed, through an analysis of the diphoton final state in data 
recorded in 2011 with the ATLAS detector at the CERN Large Hadron Collider. The analysis uses a dataset of 
2.12 fb' 1 of proton-proton collisions at *>fs = 7 TeV. The diphoton invariant mass spectrum is observed to be 
in good agreement with the expected Standard Model (SM) background. We set 95% CL lower limits on the 
scale related to virtual graviton exchange process in the context of the Arkani-Hamed, Dimopoulos, Dvali model 
(ADD) and on the lightest Kaluza Klein excitation mass in the context of the Randall-Sundrum model (RS). 



1 Introduction 

Recently, there has been great interest in models which ad- 
dress the hierarchy problem through the existence of extra 
spatial dimensions. In this analysis we search for evidence 
of extra spatial dimensions in the context of the ADD [ 1 1 
and RS |2| models. This analysis briefly summarized here, 
is described in more detail in . 

In the ADD context, n flat extra spatial dimensions are 
postulated with a compactification radius R. The gravity 
is the only field that propagates in those extra dimensions 
and acquires Kaluza Klein (KK) modes. In the ADD con- 
text, resolving the hierarchy problem implies small val- 
ues of 1/R. This leads to an almost continuous KK spec- 
trum of the graviton mass. Experimentally, an ADD sig- 
nal will contribute through virtual graviton exchange to 
the diphoton invariant mass spectrum. This process can be 
parametrized as a function of the number of extra dimen- 
sions n and an ultraviolet cutoff (M s ). Several formalisms 
exist in the literature to define M s , referred to here as GRW 
H, HLZ and Hewett @. 

The RS model postulates a 5-dimensional space-time 
bounded by two (3+1) branes with the SM particles lo- 
calized on one of the branes. The fifth dimension has a 
"warped" geometry which allows to naturally generate TeV 
scales from the Planck scale. The fundamental Planck scale 
(Mpi) on one brane is related to the apparent scale (A n ) on 
the other brane by the relation A„ — Mp\ exp (-k7rr c ) where 
k is the curvature scale of the extra dimension, r c the com- 
pactification radius and Mpi - -^=. The observed hierar- 
chy of scales can be naturally reproduced by this model 
with kr c as 1 1 . In the minimal RS model, gravitons are the 
only particles that propagate in the bulk. Consequently a 
series of massive KK excitations is predicted with a mass 
splitting of the order of 1 TeV. Finally the RS model can be 
expressed in terms of the coupling k/Mp/ and the mass of 
the lightest KK excitation (mc). Experimentally the pres- 
ence of a RS signal can be probed as a narrow resonance 
in the diphoton invariant mass spectrum. 
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2 Photon reconstruction and identification 

A complete description of the ATLAS detector can be found 
in (7). In this analysis we look for final states with two 
photons. Therefore we rely on the inner tracker and the 
calorimetric system of the detector. The tracking system of 
ATLAS is composed of layers of silicon-based and straw- 
tube detectors surrounded by a 2 T magnetic field for mo- 
mentum measurements. It allows to measure the tracks of 
the 30% of photons converting in e + <?~ pairs before reach- 
ing the calorimeter. The photon energy deposit is measured 
by the Liquid Argon (LAr) electromagnetic calorimeter, 
covering the region < \rj\ < 1.37 and 1.52 < \rj\ < 231 . 
It is segmented in three longitudinal layers. Most of the 
energy of the electromagnetic (EM) shower is recovered 
by the second layer which has a granularity of Ar\ x A<p - 
0.025 xO. 025. The first layer with a thinner segmentation is 
designed to reject jets dominated by neutral hadrons such 
as a n° or a r\. Indeed, such a hadron may decay into two 
photons which can be identified by resolving the two max- 
ima in the first layer of the LAr calorimeter. In the central 
region (I77I < 1.81), a pre sampler is used to evaluate the en- 
ergy lost in the upstream material. The hadronic calorime- 
ter, composed of scintillating tiles-iron sampling for the 
central part (|^| < 1.7) and of liquid-argon-copper/tungsten 
sampling for \rj\ > 1.7 is used to measure the leakage of 
energy of the photon candidates. A complete description 
of the reconstruction and identification of photons with 
ATLAS can be found in flS). 



3 Trigger and Data selection 

The events are recorded using a trigger requiring at least 
two photon candidates with transverse energy E y T > 20 
GeV and satisfying identification requirements based on 
the leakage of energy into the hadronic calorimeter and the 
transverse width of the EM shower. This trigger is fully 
efficient for high mass diphoton events passing the offline 
selection requirements. Events are required to have at least 
one primary collision vertex, with at least three reconstructed 
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tracks. Selected events have at least two photon candidates, 
each with E y T > 25 Ge V and lying outside the transition re- 
gion between the barrel and endcap calorimeters. In addi- 
tion, the two photons are required to satisfy standard qual- 
ity criteria and to lie outside detector regions where their 
energy is not measured in an optimal way. The two high- 
est Ej photon candidates have to satisfy a set of identi- 
fication requirements on the hadronic leakage and on the 
lateral width of the EM shower. The requirements on the 
EM shower use the thin granularity of the first sampling to 
achieve a high purity of the selected photon sample. The 
isolation transverse energy E™° for each photon is calcu- 
lated [8 1 by summing over the cells within a cone of radius 
AR = - rf) 2 + (<f>- 4> y ) 2 < 0.4 around the direction of 
the photon. Then the energy deposit of the photon itself is 
subtracted as well as the soft jet activity of the underlying 
event [9|. In addition, E"° is corrected for the leakage of 
the photon energy into the isolation ring. To further reduce 
the jet background, an isolation cut of Ej° < 5 GeV is 
applied on the two leading photons. For all events the two 
photons with the highest E y T values are considered and the 
diphoton invariant mass has to exceed 140 GeV. A total of 
6846 events are selected. 



4 Monte Carlo Simulation Studies 

Monte Carlo (MC) simulations were performed to study 
the detector response to various scenarios of ADD and 
RS models as well as to perform studies of the SM back- 
ground. 

The SM diphoton background was simulated with 
PYTHIA 6.424 (TO) and MRST2007LOMOD Hi] parton 
distribution functions (PDFs). The PYTHIA events were 
reweighted as a function of the diphoton invariant mass to 
the differential cross section predicted by the NLO calcu- 
lation of DIPHOX 1.3.2 [12 1 . The reweighting factor de- 
creased smoothly from » 1.6 for m yy = 140 GeV to 1 for 
large masses. 

SHERPA 1.2.3 [131 with CTEQ6L [ 14| PDFs was used 
to simulate the various ADD scenarios. The ADD MC sam- 
ples were used to determine the signal acceptance (A) and 
selection efficiency (e). The acceptance varied for the var- 
ious ADD implementations and fell from typical values of 
* 20% for M s = 1.5 TeV down to * 15% for M s = 3 TeV 
The selection efficiency, for events within the detector ac- 
ceptance, was found to be m 70%. 

RS model MC signal samples were produced with 
PYTHIA 6.424 for a range of mc and k/Mpi values, using 
the MRST2007LOMOD PDFs. The products of 
Axe for the RS signal models were in the range * (53 - 
60)%, slowly rising with increasing graviton mass. The 
reconstructed shape of the graviton resonance was mod- 
eled by a Breit-Wigner (BW) lineshape convoluted with a 
double-sided Crystal Ball function to describe the detector 
response. The natural width of the BW was fixed accord- 
ing to the expected theoretical value, which varies as the 
square of k/Mp/. As an example, the value of the width for 
k/Mpi = 0.1 increases from * 8 GeV for m c = 800 GeV 
to w 30 GeV for mc = 2200 GeV. Fitting this mass depen- 
dence provides a parametrization used to describe signals 
with any values of mc and k/Mp/. 



5 Background evalution 

The SM yy production is the largest background of this 
analysis. The invariant mass spectrum shape was deter- 
mined by PYTHIA reweighted by DIPHOX NLO cross- 
section predictions. 

The second significant background arises from a dif- 
ferent physics object (electron or jet) misidentified as a 
photon. This background, called reducible background, is 
dominated by the production of y + jet or dijet events. 
Studies have shown that the Drell Yan contribution can be 
neglected for events with m yy > 140 GeV. The shape of the 
reducible background is determined by a data-driven tech- 
nique. By reverting the identification criteria on the lead- 
ing and/or subleading photon, we get three different con- 
trol samples enriched in y + jet, jet + y and dijet events. 
The invariant mass spectra of the three samples are consis- 
tent and they are merged to determine the reducible back- 
ground shape. This shape is extrapolated to high diphoton 
masses where the data control sample are statistically lim- 
ited. 

The fraction of each background has been determined 
by a two dimensional template fit method using the photon 
isolation variable, described in more detail in [ 15 1. Finally 
the background expectation is normalized to the data in the 
invariant mass range [140, 400] GeV where the presence 
of any possible ADD and RS signal has been excluded by 
previous searches. 



6 Systematic uncertainties 

The systematic uncertainty on the background estimation 
arises from three different sources. The irreducible back- 
ground uncertainty is obtained by varying the scales of 
the models and the PDFs in DIPHOX while the reducible 
background uncertainty is obtained by fitting the three con- 
trol samples described in section [5] with the functionnal 
form. The uncertainty on the background estimation varies 
from x2% for the low mass region to ^20% at a mass of 
x2 TeV. For the signal, the PDFs uncertainty is of 10-15% 
for ADD models and 5-10 % for RS models and the uncer- 
tainty on the signal yields has been evaluated to be of 6.7 
% for both ADD and RS models. 



7 Results and Interpretation 

Figure Q] shows the observed invariant mass spectrum of 
diphoton events, with statistical significance of the bin-by- 
bin difference between data and the expected background 
at the bottom. This difference is measured in standard de- 
viations based on Poisson statistics. The small variations 
show the good agreement between data and the estimated 
SM background. This is confirmed by an analysis using 
the BUMPHUNTER fl]D tool, which yields a probability 
of 0.28, given the background only hypothesis, to observe 
discrepancies as large as the one observed in the bin by bin 
comparison. In the absence of any significant deviation we 
set 95% CL on the signal cross sections, using a bayesian 
approach [ 17 1 with a flat prior. In the context of ADD mod- 
els the signal search region was chosen as m yy > 1 . 1 TeV 
by optimizing the expected limit. The observed (expected) 
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Table 1. 95% CL limits on M$ (in TeV) for various implementa- 
tions of the ADD model, using both LO (k-factor =1) and NLO 
(k-factor = 1.70) theory cross section calculations. 



k-factor 
Value 


GRW 


Hewett 
Pos Neg 


HLZ 

H-3 n = 4 n = 5 n = 6 n = 7 


1 


2.73 


2.44 2.16 


3.25 2.73 2.47 2.30 2.17 


1.70 


2.97 


2.66 2.27 


3.53 2.97 2.69 2.50 2.36 



Table 2. 95% CL lower limits on mo, for various values of k/Mpi. 
The G — > 77 channel alone and the combination with the dilepton 
results of Ref. 1181 are shown, using both LO (k-factor =1) and 
NLO (k-factor = 1.75) theory cross section calculations. 



k-Factor Value 


Channel(s) Used 


95% CL Limit [TeV] 


k/Mpi Value 
0.01 0.03 0.05 0.1 


1 


G^>yy 
G -> yy/ee/iifi 


0.74 1.26 1.41 1.79 
0.76 1.32 1.47 1.90 


1.75 


G^yy 
G -> yyleelw 


0.79 1.30 1.45 1.85 
0.80 1.37 1.55 1.95 



limit on the cross section due to new physics is 2.49 (1.94) 
fb. This result can be translated into 95% CL lower limits 
on M s for different numbers of extra dimensions and for- 
malisms. Table [TJ summarizes the observed limits value. In 
the context of the RS models, the observed invariant mass 
spectrum was compared to templates of the expected back- 
ground and various signal parametrizations. The limit is set 
as 95% CL on the product of the cross-section x Branch- 
ing ratio (crB) as a function of niQ. Then the cross-section 
limit is converted into a mass limit using the theoretical 
dependence. Table [2] shows the limit in the diphoton final 
state for various values of k/Mp/. By combining with pre- 
viously published ATLAS results |18| from the dilepton 
final state, we obtain a 95% CL limit on crB as a function 
of rtiQ shown by the top plot of figure [2] The result is also 
interpreted in the plane of k/Mpi versus niQ (bottom plot 
of figure [2]). 
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Fig. 1. Observed invariant mass spectrum, superimposed with the 
predicted SM background and examples of signals for ADD and 
RS models. The bin-by-bin significance of the difference between 
data and background is shown in the lower panel. 
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Fig. 2. (Top) Expected and observed 95% CL limits from the 
combination of G — > yyjeejyiyi channels on crB as a function 
of the graviton mass. The theory bands are drawn assuming a k- 
factor of 1.75. (Bottom) The RS results interpreted in the plane 
of k/Mpi versus graviton mass, and including recent results from 
other experiments |19 20|. The region above the curve is ex- 
cluded at 95% CL. 



References 

1. N. Arkani-Hamed, S. Dimopoulos, and G.R. Dvali, 
Phys. Lett. B429, 263 (1998). 

2. L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 
(1999). 

3. ATLAS Collaboration. larXiv: 1 1 1221941 (2011); sub- 
mitted to Phys. Lett. B. 

4. G. Giudice, et al., Nucl. Phys. B544, 3 (1999). 

5. T. Han, et al., Phys. Rev. D59, 105006 (1999). 

6. J. Hewett, Phys. Rev. Lett. 82, 4765 (1999). 

7. ATLAS Collaboration, JINST 3, S08003 (2008). 

8. ATLAS Collaboration, Phys. Rev. D83, 052005 (201 1). 

9. M. Cacciari, et al., JHEP 04, 065 (2010). 

10. T. Sjostrand et al., Comput. Phys. Commun. 135, 238 
(2001). 

11. A. Sherstnev and R. S. Thome, Eur. Phys. J. C55, 553 

(2008) . 

12. T. Binoth, et al, Eur. Phys. J. C16, 311 (2000). 

13. T. Gleisberg et al., JHEP 02, 007 (2009). 

14. P. M. Nadolsky, et al, Phys. Rev. D78, 013004 (2008). 

15. ATLAS Collaboration, Phys. Rev. D85, 012003 
(2012). 

16. G. Choudalakis. larXiv:! 10 1.0390 (201 1). 

17. A. Caldwell,etal, Comput. Phys. Commun. 180, 2197 

(2009) . 

18. ATLAS Collaboration, Phys. Rev. Lett. 107, 272002 
(2011) 

19. DO Collaboration, Phys. Rev. Lett. 104, 241802 

(2010) ; CDF Collaboration, Phys. Rev. Lett. 107, 051801 

(2011) . 

20. CMS Collaboration. llarXiv: 1 1 12.06881 (201 1) (submit- 
ted to Phys. Rev. Lett.). 



